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Abstract. We search a sample of photometric luminous red galaxies (LRGs) 
measured by the Sloan Digital Sky Survey (SDSS) for a quadrupolar anisotropy 
in the primordial power spectrum, in which P(k) is an isotropic power spectrum 
P(k) multiplied by a quadrupolar modulation pattern. We first place limits on the 
5 coefficients of a general quadrupole anisotropy. We also consider axisymmetric 
quadrupoles of the form P(k) = P(fc){l + g*[(k ■ ri) 2 — ^]} where fi is the axis 
of the anisotropy. When we force the symmetry axis ri to be in the direction 
(I, b) = (94°, 26°) identified in the recent Groeneboom et al. analysis of the 
cosmic microwave background, we find g t = 0.006 ± 0.036 (Iff). With uniform 
priors on n and g* we find that —0.41 < g, < +0.38 with 95% probability, with 
the wide range due mainly to the large uncertainty of asymmetries aligned with 
the Galactic Plane. In none of these three analyses do we detect evidence for 
quadrupolar power anisotropy in large scale structure. 



1. Introduction 

Statistical isotropy (SI) is one of the most standard predictions of structure- 
formation and inflationary models. In this hypothesis, the density fluctuations in the 
Universe are a realization of a random field whose statistical properties (e.g. power 
spectra) are invariant under rotations. When probing density fluctuations using the 
cosmic microwave background (CMB) temperature, SI is generally assumed in the 
analysis. However, searches for violations of SI, or statistical anisotropy, are now 
being performed with increasing precision as the amount of CMB data has grown. 
These searches are revealing possible evidence for statistical anisotropy in the CMB 
[TJ [2] |21 S], including a possible detection of quadrupolar anisotropy [5J 0. In 
response, many have proposed inflationary and dark-energy theories with parameters 
that quantify departures from SI[S[l[Iia[IlJ[Il[ni[IS[]^ 

The detection of quadrupolar anisotropy in the CMB now appears to be contaminated 
by systematic effects (possibly beam asymmetry [23], although this is debated [B]); as 
such it is desirable to constrain quadrupolar anisotropy by other techniques. 

One way to quantify statistical anisotropy is to allow the three-dimensional power 
spectrum of matter density fluctuations -P(k) to depend on the direction of k. This 
is a full description for Gaussian but anisotropic initial perturbations. This approach 
was motivated by the inflationary model of Ref. [13], a model for which Pullen & 
Kamionkowski [24] constructed parameter estimators for CMB analysis. As usual, 
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if <5(k) is the Fourier amplitude of the fractional matter density perturbation, an 
anisotropic power spectrum is defined via 

<W(k')) = (2^) 3 <5 D (k-k')P(k), (1) 

where the angle brackets denote an average over all realizations of the random field, 
and <5rj is a Dirac (5-function; note that we still preserve the assumption that different 
Fourier modes are uncorrelated (statistical homogeneity) . A direction-dependent -P(k) 
can be decomposed via 



P(k) = P(k) 



1 + .9LAf(fc)i?LM(k) 



LM 



(2) 



where P(k) is the isotropically averaged matter power spectrum, and i?LAf(k) (with 
L > 2) are real spherical harmonics. Physically, gLM gives the magnitude of statistical 
anisotropy on order L, with M giving the direction of that anisotropy. Also, gLM = 
for odd L because a real scalar density field always has <5(k) = 6*(— k) and hence 
P(k) = P(-k). 

We define the real spherical harmonics RlmQ*) in terms of the complex spherical 
harmonics YlmQ*) by 

f ^(Ylm + Y£ m ) ifM>0 
R LM = < if M = (3) 

I lf&-^M) ifM<0; 

these are easily seen to obey the usual orthonormality rules, but have the advantage 
of making the gLM coefficients real. The expressions for L = 2 are given in Appendix 
|Appcndix A| 

The purpose of this paper is to measure or constrain the anisotropy using large 
scale structure data. Given the recent debate over the detection of quadrupolar 
anomalies in WMAP, and the evidence that the signal is contaminated by systematic 
effects [SI [7], it is worth using other datasets as well to constrain models with 
anisotropic power. In this paper we will assume for simplicity that gLM is scale- 
invariant. This is both a simplifying assumption, but is also a good first approximation 
in at least some classes of modified inflationary models [13]. We will also focus only 
on the quadrupole anisotropy g2M] this is the phenomenologically simplest type of 
anisotropy allowed, and also emerges from anisotropic inflation models in the limit of 
very weak anisotropy [15] . 

Galaxy surveys probe matter fluctuations because on large scales, the galaxy 
density is related to the matter density in accordance with a linear bias model: 

(5 g (k)d* g (k>)) = (2nfS D (k - k')^P(k) , (4) 

where <5 ff (k) is the Fourier amplitude of the fractional galaxy density perturbation, 
and b g is the linear galaxy bias. The galaxy survey probe has been used to estimate 
P(k) by stacking the measured angular matter power spectra C; in eight photometric 
rcdshift slices ranging from z = 0.2 to 0.6 [25] . By performing a similar analysis and 
including the anisotropy parameters giM in the power spectrum, we can use galaxy 
surveys to estimate quadrupole anisotropy while assuming fiducial values for the other 
cosmological parameters. 

The plan of our paper is as follows: In Section [2] we describe the SDSS data used 
and why we choose LRGs to trace the galaxy distribution. Section [3J calculates the 
angular correlations statistical anisotropy produces in galaxy surveys and constructs 
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estimators of the <72M S an d other systematic power spectrum variations. We present 
estimates of these parameters in Section^] and in Section[3]we present our conclusions. 
Wherever not explicitly mentioned, we assume a flat ACDM cosmology with 57 m = 0.3, 
Sib = 0.05, h — 0.7, n s = 1.0, and a% = 0.9. Since ours is a search for anisotropy, 
small changes in the cosmology will result only in changes in the calibration of the 
<72M estimator; they do not alter the null hypothesis. 

2. Choice of sample 

There are several ways to use galaxy survey data to search for statistical anisotropy. 
In principle, one could use a 3-dimensional redshift survey and search for anisotropy 
in the power spectrum. This would however be very technically involved: redshift- 
space distortions make the line of sight direction special. With sufficient sky coverage 
one could break the distinction between redshift-space distortions and true statistical 
anisotropy. However, in this paper we choose the technically simpler route of using 
photometric galaxy catalogues, which can be studied using estimators analogous to 
those for the CMB. 

The photometric data we use come from the Sloan Digital Sky Survey (SDSS) [26 . 
The SDSS consists of a 2.5 m telescope [27] with a 5-filter (ugriz) imaging camera [28] 
and a spectrograph. Automated pipelines are responsible for the astrometric solution 
[29] and photometric calibration (30j [31j [32j [33] - Bright galaxies, luminous red galaxies 
(LRGs), and quasars are selected for follow-up spectroscopy [34j [35] [36] [37] . The data 
used here were acquired between August 1998 and October 2004 and are included in 
SDSS Data Release 5 [38]. 

We use a sample of photometrically classified luminous red galaxies. LRGs 
are the most luminous galaxies in the universe, making them appealing for probing 
maximal volume. They also tend to be old stellar systems with uniform spectral energy 

distributions and a strong discontinuity at 4000 A, which enables precise photometric 
redshifts and hence measurements of {g2ki} in multiple redshift slices. This both 
reduces statistical error bars and allows tests for consistency. The cuts that define 
the photometric LRG sample are enumerated in Ref. [39J. The sample is divided into 
8 photometric redshift slices of thickness Az p = 0.05 ranging from z P! min = 0.2 to 
%>,max = 0.6, using the "single-template" photo- z algorithm of Ref. [39]. We plot the 
redshift distributions in Fig. [TJ while their properties are given in Table [1] 

Our galaxy catalogue is a subset of that used by Ref. [ID], but we restrict our 
attention to Galactic latitudes b > 45°. This was done in Ref. 25 to prevent stellar 
contamination in the data, and we decided to use the same cut accordingly. We pixelize 
these galaxies as a number overdensity, S g — (n — n) /n, onto a HEALPix pixelization 
[4T] of the sphere, with 1,418,213 pixels. This corresponds to a solid angle of 4662 
deg 2 (as opposed to 3528 used in Ref. [H]). The LRG maps thus generated are shown 
in Fig. H 

3. Formalism and estimators 

3.1. Galaxy density projections on the sky 

We relate angular correlations in the sky to the direction-dependent matter power 
spectrum -P(k), which follows the derivation in Ref. [25] for the statistically isotropic 
case. A photometric galaxy survey measures the 2-dimensional projected galaxy 



Anisotropic power spectrum 4 



20 




z 



Figure 1. The rcdshift distributions for the 8 photometric redshift slices. 




Figure 2. The LRG density in the 8 photometric redshift slices. The 45° radius 
caps are displayed in a Lambert Azimuthal Equal- Area Projection, with the North 
Galactic Pole at the centre, I = 0° at right, and £ = 90° at bottom. The labels 
indicate the characteristic redshift of each slice. 



Anisotropic power spectrum 



5 



Table 1. Properties of the 8 redshift slices; z p is the photometric redshift range, 
and z meaI1 is the mean (true) redshift of the slice. N ga i is the number of galaxies 
in the redshift slice, and b g is the linear bias. 



Label 


Zp 


2 mcan 




bg 


zOO 


0.20-0.25 


0.233 


30812 


1.74 


zOl 


0.25-0.30 


0.276 


38168 


1.52 


z02 


0.30-0.35 


0.326 


37963 


1.67 


z03 


0.35-0.40 


0.376 


55951 


1.94 


z04 


0.40-0.45 


0.445 


77798 


1.75 


z05 


0.45-0.50 


0.506 


138901 


1.73 


z06 


0.50-0.55 


0.552 


126318 


1.80 


z07 


0.55-0.60 


0.602 


93973 


1.85 



overdensity <5 g (n), which is related to the full 3-dimensional density via 

S g (n) = J d X /(x)5 g (x = X ft), (5) 

where x 1S the comoving distance, and /(x) is the radial selection function, which is 
the normalized redshift distribution for the redshift slice. (We leave out redshift space 
distortions here, but include them below.) By Fourier transforming S g (x) and using 
the identity 



e -ikflx 



^(2;+l)i^(fcx)fl(k-fi), (6) 

1=0 

we obtain 

oo /* j3 7 

5 3 (fi) = ^i i (2/ + l) J — i^n^k)^*), (7) 

where 

Wi{k) = J d x f(x).n(kx) (8) 

is the window function, and j\ (x) and Pi (x) are the Z^ n -order spherical bessel functions 
and Legendre polynomials, respectively. 

The statistical properties of the 2-dimensional galaxy field can be derived in 
analogy to those for the CMB [21], with Wi(k) replacing the CMB radiation transfer 
function If statistical isotropy is valid, then the two-point galaxy correlation 

function can be written as 

C g (h, n')|si = (S g (A)5 g (h')) |si = £ — — C 9 , ; P ; (n • ft') , (9) 



47T 
i 



where C g ^i is the angular galaxy power spectrum: 



C g ,i = - I dkk 2 P g (k)[W,(k)} 2 , (10) 



7T 



I Note that with our definition of the Fourier transform, there is a relative factor of (27r) 3 between 
some of our formulas and those of Ref . |24| . 
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where P g (k) = b 2 g P(k). However, the presence of statistical anisotropy will cause 
additional terms to appear other than the ones in Eq. |H1 Using Eqs. ([7]) and @, we 
find for the statistically anisotropic case 

i 

+ E E D^x^ m 'Rumi^(&'). (ii) 

LM Iml'm' 

Here, the set of C g ^s are given by the usual galaxy power spectrum for the case of 
statistical isotropy. Statistical anisotropy produces the second term, where 

/ 2 f°° 

D g ,u'=i l ~ l - dkk 2 P g {k)g LM Wi{k)W v {k), (12) 
77 Jo 

and 

Xlml'm>= J #&Rlm$)Rv m '$)RLM$) (13) 

is the real 3-harmonic coupling coefficient. Its expression in terms of Wigner 3-j 
symbols is given in | Appendix B| Parity implies that X^f, m , is nonzero only for 
I + I' — L even. Since earlier in the Introduction we showed L is even, this requires 
l + V to be even and D^^, to be real. Equations (fT0|) and (fl2j) agree with similar results 
in Ref. [23], and X^,' m , in Eq. (H3J is analogous to in Ref. [23 ]. Throughout 

this paper, we use upper-case indices LM for statistical anisotropies in the matter 
power spectrum, and lower-case indices Im for random anisotropies in the galaxy 
distribution. 

In Fig. |3jwe show the predicted angular galaxy power spectra C g j for the eight 
redshift slices in our analysis assuming our fiducial cosmology. We use CMBFast [32] 
to calculate P(k), and we use the Halofit prescription [33] to evolve P(k) into the 
nonlinear regime. Note that we do not attempt to account for the nonlinear evolution 
of QlMi which only suppresses the primordial anisotropy by < 7% on quasilinear 
scales [44]. We also display [— F g ^i +2 )], where F gt u> = D g ^l / Qlm- We only need 
to show F g j(i+2) because we are interested only in quadrupolar statistical anisotropy, 
or L = 2. For this case it can be shown that Xf^j,^ is zero except for the cases 
V = I and V = I ± 2, and F g ^i = C g .;I§| Notice that for large I and smooth f(x)i 
Wi{k) -> ^Jlf(l/k)/(2k); then we have W l+2 ~ Wi and F gMl+2) ~ -C g , h 

We also show in Fig. [3] C g j and [— F g i^ +2 )] when the effect of redshift space 
distortions is included. We include this effect by substituting for the window function 
Wi = + Wf, where W® is the window function shown in Eq. [S]and W[ is given by 



W[(k) = (3 



2/ 2 4 9/ — 1 



(2/ + 3)(2Z-l) 
1(1-1) 



;Wl 2 {k) 



(2/-l)(2/ + l) 

(14) 



(2Z + l)(2Z + 3) Wl+2{ ' 



§ For nonzero X^J^, ,, the triangle inequality requires \l — l'\ < L = 2, and parity requires I — l' to 
be even. 
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and j3 is the redshift distortion parameter given approximately by j3 ~ fl^ /b g . The 
formulas for redshift space distortions in the angular galaxy power spectrum were 
derived in Ref. [25] . 

3.2. Estimation of power and statistical anisotropy 

We construct a quadratic estimator HB] to measure the anisotropy coefficients 
92M- As always when searching for anisotropics, it is necessary to fit simultaneously 
for the galaxy power spectrum, the anisotropy, and any systematics terms that may 
be present in the data. The basic premise is to write the galaxy density fluctuation 
map as a vector S g of length N p - lx . This vector has an iV p i x x N p i x covariance matrix 
C, which we parameterize as 

N t 

C = ^ Pi C i; (15) 

i=l 

where the {pi}^J; 1 are parameters to be estimated and {C l j}^ 1 are "templates." The 
notation also serves to remind us that = DC/dpi. In the case of a Gaussian 
random field, the minimum- variance unbiased quadratic estimators for the {pi} are 

p = F- 1 q, (16) 

where 

F« = itr[C, i wC, j w] (17) 

is the Fisher matrix and 

ffi = i^wC,(W<5 g . (18) 

Here w is a weighting matrix, which should be taken equal to the inverse of the 
covariance matrix C. Fortunately, the estimator p remains unbiased (but not 
necessarily minimum-variance) for any choice of weight w, and regardless of whether 
the true galaxy field 6 g is Gaussian or not. For our analysis, we take the weight to be 
w = (S + N) _1 , where S is the signal covariance matrix (diagonal in Zm-space, and 
using the theoretical power spectra) and N is the Poisson noise. The matrix inversion 
and trace estimation are done by the iterative and stochastic methods described in 
detail in Refs. g2|25]. 

We next turn our attention to the template construction. The simplest template 
is the Poisson noise itself, 

Cy.iv = "_r , (19) 
n 

where n is the mean number of galaxies per pixel, and N is the noise amplitude (1 
for pure Poisson noise). We may also parameterize the isotropic part of the power 
spectrum by band power amplitudes C n with C g j = Yl n =i Cnifi, where 77" is a step 
function that is 1 when I is in bin n and otherwise. The corresponding template is 

-Jl=Y / Rlm(n i )Rl m (n j )r 1 r- (20) 
!b Im 

We use 18 bins in I, ranging from I = 2 up through 600. 

For the anisotropy parameters giu-, the templates are somewhat more 
complicated. We first extract g2M from the anisotropy amplitude D 2Ai ■ 



g.lV 

Dg,U> = 92MFg.IV , (21) 
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Figure 3. The predicted angular power spectra for each of the 8 redshift slices. 
The solid lines show nonlinear auto power spectra, while the dash-dotted lines 
show the cross power power spectra with the adjacent slice at higher redshift. 
The dashed lines show the predicted [—F g< i(i+2)]- The dotted lines show the 
effect of redshift space distortions on all three spectra. 
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F, 



g,W 



■l-V 



&kk 2 P g (k)Wi(k)W v {k). 



(22) 



Inspection of Eg. [TTIthen leads us to: 



dC, 



dg 



'2M 



E 

Irara' _ 



F g,l{l+2) X hn\l+2)rn' R lrn{^i)R(l+2) m '{^j) 
Fg,l{l-2) X lm(l-2) m > R lrn (fij )R(l-2)m> ) 



(23) 



In addition to these "essential" templates, we have also added two others to 
project out various systematics that could mimick statistical anisotropy. In particular, 
if there are photometric calibration errors that vary slowly across the survey (either 
coloured or grey then the depth or effective redshift may vary, which would lead to 
a modulation of both the signal power spectrum and the noise level across the sky. 

We model the modulation of the signal power spectrum by considering a 
modulation in the fractional density perturbation field in the form <5'(n) = [1 + 
h(h)]S(n). This modulation will cause the two-point galaxy correlation function to 
have an extra factor of [l + /i(n)][l + /i(n')] ~ l + h(h) + h(h'), assuming h(h) <C 1. We 
choose to allow h(n) to have a quadrupole pattern. (Other forms of slow modulation 
across the sky, e.g. a dipole, should be degenerate with a quadrupole given that our 
data is only in a cap of radius 45°.) 

In this case S'(h) can be written as 



S'(h) 



1 + 



2 

E 

M=-2 



^2A/-fiW(n) 



5(h) 



(24) 



where hiM are the modulation parameters. In the case of modulation with no 



statistical anisotropy, C 3 (n, n') is given by 



C 9 (n,n')= 1 



(n, n')|s/ , (25) 

M=-2 ) 



where C 9 (n, n')\si is given by Eq. ©. By using the identity 



Rl m (h)Rl> m '{h) - ^ ^Irfd'-m' Rlm{^) , 
LM 



(26) 



we find 



C 9 (n,n / ) = C 9 (n,n / )|si 

h 2M (C g ,i + C g , v )X^ Vm ,Ri m {h)R Vm ,{h'),{27) 



H E 

Iml'm' M 



hence 



dh 



2M 



E 

Imm' 



2Cg,l X lmlm' Rimini) Rim' 



|| Coloured errors apply to an error in the relative calibration of different bands, e.g. g — r, whereas 
grey errors leave colours fixed but vary the magnitude of an object. 



Anisotropic power spectrum 



10 



+ (C g ,i + C a ^ +2 )Xf^ (l+2)m ,Ri m {^i)R(i+2) m '{^]) 

+ (Cg.l + C s ,;- 2 )^ 2 r ^(i_2) m '^m(n0^(i-2)m'(nj) 

An analogous construction for modulation of the Poisson noise gives 



df 2 M n 

These 10 parameters (h 2 M and J 2 m) are jointly estimated with {C n , <72m}- 



(28) 



(29) 



3.3. Gaussian simulations 

We test our estimator on a suite of simple simulations in order to verify its ability 
to detect anisotropy when it is present (and to measure zero when anisotropy is 
not present). Gaussian simulations are sufficient for this purpose since a quadratic 
estimator by construction cannot be sensitive to higher moments of the data. 

We perform two tests, one without anisotropy or modulation and one with both. 
In each test, we use the power spectrum C g j and F g m+ 2 ) of redshift slice zOO for our 
fiducial cosmology to construct 50 sets of simulated galaxy perturbation maps over 
the pixels in our analysis' viewing region. We also add Gaussian noise to each pixel, 
with the variance in the noise set equal to 1 /n ga i , where n ga i is the average number of 
galaxies per pixel for the redshift slice (consistent with Poisson fluctuations). Then, 
we run each simulation through the algorithm to find an estimated set of parameters 
C n , g 2 M, and h 2 M- We then average these parameters over all 50 simulations to find 
an output set of parameters, which we then compare to our input parameters for 
constructing the simulations. Since we do not input C n directly, we instead compare 
the output C n to Cg^i at the median I of bin n. Note that the variance used to compare 
the input and output parameter sets is equal to the variance of one simulation, taken 
from the diagonal of the inverse fisher matrix, divided by the number of simulations. 

For the first test, our simulations had input parameters g 2 M — h 2 M — for all M . 
Since there is no covariance between the simulated 5 g j m , the real spherical harmonic 
coefficients of S g (n), we can simulate each 5 g .i m independently. A plot of the input 
and output parameter values for C n , g 2 M, and H 2 m are shown in Fig. [4] In the figures 
we see good agreement between the input and output values. This test shows us that 
our algorithm should not see anisotropy or modulation where there is none. We also 
see that the error for c/2,0 is larger than the errors for the other g 2 MS. This is due to 
a lack of data in the equatorial plane while (72,0 parametrizes a quadrupole along the 
z-direction. Our data set makes this type of quadrupole less distinguishable from a 
uniform excess over the whole sky than with other quadrupole types. 

In the second test, we simulate anisotropic power and modulation in the z- 
direction by setting 1720 = ^20 = 0.5. Now that the <5 ffj ; m s are correlated, their 
simulation is no longer trivial. To construct the simulation, we define the matrix 
Ej;™\ equal to (5 g j m 6g,i'm) in the case g 2 u — h 2 M — f° r all M except M — 0. This 
matrix is given by 

E i™ ) = CgjSw + g 20 F g:W Xf° nl , m + h w (C gt i + C gt i>)X?° d , m . (30) 

To construct our simulation, we perform a Cholesky decomposition on E^ m ^ to find 
the triangular matrix L< m ) such that E^™) = L ( - m ^L^ m '> T . We use L/ m ) to construct 
S^"^ = l/ m )x( m - ) , where x( m ) is a Gaussian random matrix with zero mean and unit 
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J 0.2- - 



-0.2 I , , , i , , , i , , , i , , , I 

-4 -2 2 4 

M 

Figure 4. The parameter values for the simulation test with no input anisotropy 
or modulation, including 1-sigma errors. The top left panel shows the C n , the 
top right panel shows the <72M s i an d the bottom panel shows the /i2M$- 

variance. This setup gives us the desired covariances we need in the S g .i m s. A plot of 
the input and output parameter values for C n , 92M, and hiM in this case are shown 
in Fig. [SJ In the figures we see good agreement between the input and output values, 
including for (720 and ^-20- 

4. Results 

We show our results for C g j in Fig. [5] Because our galaxy sample covers only 11% 
of the sky, the powers in C; and giM become degenerate, and the isotropic and 
anisotropic templates become highly correlated. Thus, when calculating C;, we remove 
the anisotropic templates for g 2 M, h 2 Mi and / 2 a/ so that power in the C;S is favored. 
When using this approach, we see agreement between the measured power spectra 
and the predicted spectra. 

Our results for giu are shown in Fig. [7J For each multipole of g 2 M, we see 
consistency with the null result among the redshift slices except for the measurement 
of 320 m redshift slice z03. At this redshift, we measure 320 = 0.925 ± 0.258 (Fisher 
uncertainty) or ±0.315 (uncertainty derived from iV-body simulations, as described 
in Sec. 14. ip . This formally corresponds to a 3. 59c (Fisher) or 2.94c (simulation) 
detection significance; however all of the other redshift slices have 320 within lcr of 
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Figure 5. The parameter values for the simulation test with input anisotropy 
and modulation, including 1-sigma errors. The top left panel shows the C n , the 
top right panel shows the </2M s > and the bottom panel shows the /l2M s - Note 
that the input values for 320 anQl ^20 are 0.5. 



zero. This is puzzling and in principle could indicate either a statistical fluke or a 
systematic error that afflicts the z03 slice. We note that the statistical significance 
is marginal: given that we calculated 5 x 8 = 40 (72A/S, the probability of having at 
least one of them deviate by 2.94cr is 12% (assuming a Gaussian distribution). On the 
other hand, the z03 slice is also the redshift at which the LRG colour locus changes 
direction [25]. The z03 slice also has the highest bias, which would make is susceptible 
to nonlinear errors. 

The results we found for the other multipoles were consistent for each redshift 
slice only when we allowed Kim and Jim to vary from the null result. We show the 
results for h^u m Fig. [8] and f^M m Fig- El Note that many of the /12MS and /2MS 
are inconsistent with zero, which hint at possible systematic errors of these forms. 



4-1. Combined statistical anisotropy estimate 

To find an estimate of g^M combining all of the redshift slices, we construct a minimum- 
variance estimator of the form 



92 M 



T,i92M,t/0- 



92M.. 



2 

92 M,- 



(31) 




Figure 6. The measured angular power spectrum for the eight redshift slices. The 
solid lines are the predicted nonlinear power spectra for our fiducial cosmological 
model, and the crosses are the measured spectra. 
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Figure 7. The quadrupole anisotropy parameters vs. rcdshift slice for each 
multipole with lcr errors from the simulations. Note 320 = 0.925 for rcdshift 
slice z03, formally a 2.94<r detection. 



where <?2A/.i is the estimate of g2M in redshift slice i and a g ^ u i — (F )mm for redshift 
slice i. A crude estimate of the uncertainty in g~2M is given by 

^ = S^- < 32 > 

92M i Q2M,i 

This uncertainty estimate neglects covariances between the redshift slices and non- 
Gaussian (trispectrum) corrections to the errors in individual slices. We therefore 
expect it to somewhat underestimate the true uncertainty in <?2M ■ For this reason we 
expect that error bars derived from iV-body simulations (as described next) are more 
reliable. The g^M values and their uncertainties as calculated by Eq. (|32|) are shown 
in the top panel of Fig. [101 

We may alternatively estimate the covariance matrix Cmm' of g2M using TV- 
body mock catalogues, which contain the correct slice-to-slice correlations and a more 
realistic description of the true non-Gaussian density field. We used a suite of 10 
simulation boxes of size (2/i _1 Gpc) 3 with periodic boundary conditions, described in 
more detail in Ref. |48j . For simplicity, and since our objective is to obtain a covariance 
matrix rather than a precision prediction of the power spectrum, we have used the 
halo catalogue from a single simulation output at z = 0.3. We populate each halo with 
a galaxy (or two galaxies if Mhaio > M 2 ) and use its "true" redshift (including the 
halo peculiar velocity) and the photo-z error distribution [39] to assign a photometric 
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Figure 8. The quadrupole modulation parameters vs. redshift slice for each 
multipole with la errors from the simulations. 



redshift. A catalogue of galaxies is then constructed by taking each halo down to some 
minimum mass M min j fixed by the requirement to have the correct number of galaxies 
in the ith photo-z slice. The parameter M2 controls the amplitude of the "1-halo" 
term in the power spectrum arising from multiple galaxies per halo (in the sense that 
the 1-halo term is set to zero if Mi = 00). We choose M2 by first constructing a mock 
catalogue with M2 = 00. The excess power AC; in the 300 < I < 600 range is then 
determined for each photo-z slice. The 3D number density of doubly-occupied haloes 
rid required to produce this excess power is then estimated as 

fi 2 

nd = 2dv/dn ACl - (33) 

An average value of rid is taken over all slices (rid = 1-8 x 10 _5 /i 3 Mpc -3 ) and this is 
used to set a mass threshold (M 2 = 7.8 x lO 13 ^" 1 A/ Q )§ This procedure generates 
an entire simulated photo-z survey, including all 8 slices and the correct correlations 
among different slices due to large scale structure. 

We construct 40 realizations of the survey, using each of the 10 boxes 4 times with 
different observer locations. We then estimate g^M for each simulation and redshift 
slice, which we marginalize using Eq.|3T]to find an estimate of t/^M f° r eacn simulation 



K In principle, this procedure could be iterated with computation of a new AC;, etc., however the 
method is probably too crude to justify such a procedure. 
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Figure 9. The Poisson noise modulation parameters vs. rcdshift slice for each 
multipole with lc errors from the simulations. 



Table 2. The anisotropy coefficients, averaged over the 8 redshift slices, using 
the more conservative covariance matrix from the TV-body simulations. Their 
covariance matrix is also given. 



M 


92A//10" 2 




Cov[32M,92M']/10 3 




-2 


3.90f 


1.238 


-0.327 


0.741 


0.879 


-0.283 


-1 


-0.979 


-0.327 


3.599 


-0.164 


-0.566 


-0.042 





9.508 


0.741 


-0.164 


11.866 


-0.761 


-1.348 


1 


6.479 


0.879 


-0.566 


-0.761 


5.044 


-0.015 


2 


-2.235 


-0.283 


-0.042 


-1.348 


-0.015 


1.106 



a. We use these estimates to construct the covariance matrix of the form 

1 40 

C M M> = 7^ (#2M - 92M)(92M' ~ 5 2 M' ) . (34) 
a— 1 

where g 2M is g 2 M averaged over the simulations. The diagonal elements of Cmm> give 
the uncertainties in 32 m- We plot 172 m with these uncertainties in the bottom panel 
of Fig. [TOl We see in this case all the measurements are within one sigma of the null 
result, which is consistent with statistical isotropy. 

The final values of giu and their covariance matrix are given in Table [5J 
Finally, as a systematics test, we consider how the <?2A/s change when we do not 
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Figure 10. The quadrupole anisotropy parameters for each multipolc 
marginalized over redshift slice with 1-sigma errors. The top panel includes errors 
calculated from the Fisher matrix. The bottom panel includes errors calculated 
using N-body simulations. Note that both results are consistent within two sigma 
with the null result, shown as the dashed line in both plots. 



project out the extra templates {/2M, ^m}- We have shown in Table[3]the changes in 
92M when none of the systematics templates are included ( "<72M only" ) and when the 
/ 2 m templates are left out but h 2 M is included ("<72M & ^2m")- As we can see from 
the table, the exclusion of the / 2 m templates has essentially no effect, but there is a 
substantial change in 52M when the h,2M templates are excluded as well. However, 
since we expect the main effect of systematic power spectrum modulation across the 
sky to be taken into account via the /12A/S, and given that they change the result by 
< 3cr, we do not expect a significant residual systematic after the /12A/S and /2MS have 
been projected out. 

4- 2. Comparison with CMB results 

Groeneboom et al. [6] report evidence for a quadrupolar power asymmetry in the 

5- year WMAP data. They investigated models of the form: 

P(k)=P(fc)[l+.g»(k-n) 2 ], (35) 
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Table 3. The changes in the anisotropy coefficients, averaged over the 8 redshift 
slices, when only some of the systematics templates are included. We present 
A;?2Af values, which are equal to the j2M from the full analysis minus those with 
only some of the systematics templates, and the number of sigmas by which the 
correction differs from zero, AgsM /^{^92m). 







92M only 




92M & h,2M 


M 






(S2M) 




-2 


0.0128 




0.63 


0.0004 


-1 


0.1591 




2.02 


0.0027 





0.2827 




1.94 


-0.0050 


1 


0.2830 




2.78 


0.0004 


2 


-0.0003 




0.02 


-0.0003 



where g* is the amplitude of the asymmetry and n is its preferred axis. For <C 1, 
this is equivalent to our Eq. ([2]) with 



LM 



9lmRlm{^) = g* 



(k-fi) 2 --; 



^P 2 (k-n) (36) 



and a slightly rescaled definition of the power spectrum, P(k) = P(k)(l — |<?*)- Here 
the — ^ ensures that there is no L = term in Eq. (j36]) . Using the spherical harmonic 
addition theorem, we can see that this requires 
Stt 

gLM = —rg*S L 2R2M('n). (37) 
15 

We may use this to construct an estimator for assuming a particular direction n; 
here we will take n to be in the Groeneboom & Eriksen direction so that we can test 
for consistency with their value of <?* . The best estimator is 

- _ 15 Emm' [^°' v ~ 1 ]mm'92mR2M' (n) 

MM' [Cov ]MM'-R2Ai" (n)i?2M'(n) 

with uncertainty 

a&) = ^.-— 1 ; (39) 

yJ2MM'l C ° V 1 }MM'R2M{n)R 2 M'(fl) 

here Cov is the 5x5 covariance matrix of the estimators for g2M (see Table [2]). 

Using the WMAP W-band maps, and considering multipoles in the CMB out to 
; max = 400, Groeneboom et al. [5] find an asymmetry of 3^ MB = 0.29 ± 0.031 with 
an axis of maximum power in the direction (l,b) = (94°, 26°), which they attribute to 
an unknown systematic effect because different signals are observed in the V and W 
bands and the apparent alignment with the Ecliptic Poles. Using the above projection 
procedure, we find an amplitude g* RG = 0.006 ± 0.036 in this direction. Groeneboom 
et al. also did their fit using the WMAP V-band maps, finding g™ B = 0.14 ± 0.034 
in the direction (l,b) = (97°, 27°); when we project our LRG anisotropy coefficients 
onto this axis, we find g^ RG = 0.007 ± 0.037. Foregrounds and noise mis-estimation 
have been disfavored as possible candidates [7] . A possible cause for the appearance of 
statistical anisotropy in the CMB data would be the ellipticity of the WMAP beams, 
which when combined with the survey strategy could result in a preferred axis in the 
direction of the Ecliptic Poles [5J HH1 US! 123] ■ Specifically, Hanson et al. [23 find that 
once asymmetric beam effects are subtracted, the data is consistent with the isotropic 
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model; however Groeneboom et al. [5] evaluated the resulting effect and found it to 
be negligible. The WMAP 7- year analysis finds no known instrumental effect other 
than beam asymmetry that can cause the anomaly, but they have not yet completed a 
full simulation of beam asymmetry effects on quadrupolar power modulation [7] . We 
also note the the WMAP team has already accounted for these beam effects in their 
estimation of the power spectrum, so this systematic in the quadrupolar anisotropy 
does not affect the cosmological parameters derived from WMAP. 

Thus the cause of the apparent asymmetry in the WMAP maps is not definitively 
known. In any case, our LRG analysis finds no anisotropy in this direction. It is 
possible that g 2 u is different at the two scales probed by the CMB and the LRG 
sample. In most variants of inflation, where each e-fold of expansion is similar to the 
previous one with ~ 0(1/ N) deviations (where N is the number of e- folds remaining), 
we would expect g 2 M to vary smoothly with the number of e-folds, or In k. We find the 
effective scales probed by Groeneboom et al.'s CMB analysis [5] and our LRG analysis 
are 0.020 Mpc" 1 and 0.15 Mpc -1 , respectively (see | Appendix C ), which differ by only 



2.0 e-folds. It would be very surprising if inflation were not only anisotropic but also 
managed to produce a scale-dependent anisotropy that varied over so short a baseline. 

4-3. Direction-marginalized constraint on 

The above analyses have either set constraints on a general g 2 u (a 5-dimensional 
parameter space) or on g* for a fixed anisotropy axis (a 1-dimensional parameter 
space) . It is however of interest to set constraints on general axisymmetric quadrupolar 
anisotropies, such as Eq. (|36|) . which would arise if there were a single preferred axis 
during inflation. This is a 3-dimensional parameter space: there is an amplitude g* 
and a direction n eS 2 . 

We may set constraints on <?* via a Bayesian analysis in which a uniform prior is 
placed on n, as has been done in several previous statistical anisotropy analyses [5ll51j. 
Our problem - setting a limit on the amplitude of an anisotropy while marginalizing 
over its direction - is similar to that performed by Ref. [51| for the large scale structure 
dipole; we follow the same methodology, although we note that for the quadrupolar 
asymmetry g* could be positive or negative ( "prolate" and "oblate" power anisotropies 
are different and cannot be rotated into each other). The marginalized likelihood 
function for is 



= /expi--^[Cov ^mm- 

J I MM' 



92M - —-g*-ti 2 MW 

15 



92M> - —g*ti 2 M>{n) 
15 



d 2 n, (40) 



where g 2 M are the estimated anisotropy coefficients and Cov is their 5x5 covariance 
matrix. If we set a uniform prior on g*, as done by Groeneboom & Eriksen [5], then 
we may divide £(<?*) by its integral J C(g*) dg* and treat it as a posterior probability 
distribution. If we do this, then we find that 68% of the posterior distribution is 
contained within -0.12 < g* < +0.10 and 95% within -0.41 < g* < +0.38. Note that 
the distribution has very non-Gaussian tails because of the large uncertainty on g 2 o : 
a quadrupole anisotropy aligned with the Galactic axis would be difficult to detect 
given our sky coverage. There is a small probability for such an alignment to occur 
and not produce measurable g 2 u (M ^ 0) even if g* is large. 
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5. Conclusions 

We have conducted a search for statistical anisotropy in the galaxy distribution. 
Statistical anisotropy can manifest from the direction-dependent primordial power 
spectrum shown in Eq. [5] with the magnitude of the anisotropy parametrized by gLM- 
This phenomenon causes the angular galaxy power spectrum C 9j ; to be generalized by 
Dgjf,, which includes gLM- We used estimators formulated by Padmanabhan et al. [25] 
and a sample of LRGs from SDSS to search for evidence of quadrupolar anisotropy 
parametrized by giM- We found g2M for all M to be within 2a of zero. Using our 
estimates of g^M and assuming a symmetry axis in the direction (l,b) = (94°, 26°), 
we calculated the anisotropy amplitude g* = 0.006 ± 0.036 (la). This confirms that 
the previously identified anisotropy in the WMAP maps (already believed to be a 
systematic effect) is not of primordial origin. When marginalizing over the symmetry 
axis direction and assuming a uniform prior for g*, we constrain —0.41 < g* < +0.38 
with a 95% confidence level. 

Looking forward, we expect much better sensitivity to from future galaxy 
surveys. For fixed sky coverage, the uncertainty in giM is proportional to the inverse 
square-root of the number of modes measured, i.e. it is proportional to l^xN z 
where /max is the maximum multipole at which the galaxy distribution is well-sampled, 
and N z is the number of effectively independent redshift slices. The largest advance 
may be possible with future large-volume spectroscopic surveys intended to study 
baryon oscillations. Here the effective number of redshift slices is N z ~ /c max Ar/7r, 
where Ar is the radial width of the survey; for surveys that reach out to z f=s 2 this is 
N z ~ 100 (instead of 8 here). As this redshift corresponds to a factor of ~ 3 increase in 
distance relative to the SDSS LRGs, we would expect that for similar sampling nP(k) 
^max should increase by a factor of 3. Thus such a survey should in principle be able 
to improve measurements of giu by an order of magnitude relative to those presented 
here. Further improvements in (720 might also be possible if improvements in the dust 
map or work in redder bands allows one to work at lower Galactic latitudes. 
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Appendix A. Real / = 2 spherical harmonics 

In Eq. [3] we introduce our convention for the real spherical harmonics Rlm (&,<]>)■ To 
clarify the functional form of Rlm, wc list the harmonics for L = 2. These are given 
by 

ifcaCM) = sin 2 9 cos(20) 

R-)\{9,4>) = —\ — cos 9 sin 9 cos <j> 

V 47T 



R w {9,4>) = ^jj- (3 cos 2 9-1) 
R2 -i(9, 4>) = \ — cos 9 sin 9 sin (f, 

V 47T 



R 2 ,- 2 (9,4>)= - W-^-sin 2 0sin(20). (A.l) 

V 1D7T 

Appendix B. Expressions for the anisotropy coefficient 

In Ref. [24], Pullen and Kamionkowski introduced an anisotropy coefficient Ciw^m' 
that appears in the correlation function, given by 



& m - = J dky; m (k)^, ra ,(k)7 iM (k) 



= {-l) m (Gt v ) 1/2 C^,._ m ,, (B.l) 
where Cf^, m , are Clebsch-Gordan coefficients, and 
L _ (2/ + l)(2Z' + l) , L0 > a 
G "' = 4tt(2L+1) [Cm ' o) ■ {B - Z) 
However, since we use real spherical harmonics (given by Eq. [3]) in our analysis 
as opposed to complex spherical harmonics, we introduce the anisotropy coefficient 
-^imi'm' gi yen by Eq. (TIB")) . We choose to write X^, m , in terms of Wigner 3j symbols. 
Due to the piecewise nature of the real spherical harmonics, X^, rn , will have different 
expressions for different values of m, m', and M. After much algebra, we can find the 
expressions for X^f, , in terms of Wigner 3j symbols (written in matrix form) and 
Pwl, given by 



/ (2i + l)(2Z' + l)(2L + l) ( I V L , 
Pwl = \I ^ I' (R3) 
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which is nonzero only for I + I' + L even. The expression for M = is given by 

XZ'm' = {-l) m PwL^ l m V m q)w. (B.4) 
The expressions for M ^ can be obtained similarly, e.g. 

~(-l) m ' (I V L 



-vLM p 



^2 V m — (fn + M) M 
(-l) m / i V L 



8m' 



m' ,m+ M 



Jm'm—M 



(B.5) 



\ "i M — m —M 
for m > M > 0; the other equations are similar but will be omitted for brevity. 

Appendix C. Effective scale for quadrupole asymmetry analyses 

In Ref. [BJ, Groeneboom et. a!, calculated a quadrupole asymmetry in the matter power 
spectrum by analysis of the CMB up to multipoles of I — 400. In this appendix, we 
derive the effective wavenumber k e g of this CMB measurement, as well as fc e ff for our 
measurement using the LRG sample. 

To find fc e ff for the CMB analysis, we first find an estimator for g 2 M in terms of 
measurable quantities in Fourier space0 This calculation has been done previously 
in Ref. [24] . By using Eq. (23) in that paper, we can construct the minimum- variance 
estimator for g%M given by 

a _ Zi,hM,u>Qu'2(F™n 2 /(cr p crn 

where C ; map = | | 2 C; + C" is the map's power spectrum, Wj = e~' 2<T ^ 2 is the 
beam window function, F™ p = W x W v F lv , Q W2 = Ew Www) 2 = (fli'2) 2 /5> and 
g2M,w is the estimator of g2M- We can construct an estimator for each IV pair with 
g2M,w = DfjP/Fui, where Dfft 1 is an estimator constructed from the measured ai m s. 
By approximating for large I 

I I 2 \ 2 1 , / I l±2 2 \ 2 3 



J ~ 8/ { ) ~ 167 ' (C ' 2) 

using Eq. (8.32) of Ref. [52], along with F u = Q and F l(l±2) ~ -0.5C, (for 
temperature perturbations), we have 

Y<1 lTlls^92M,ll + 64^(?2M,i(i+2) + ff2M,I(I-2))] , „s 

52M ^ 7 , 

32^1^l ll l 

where T t = [1 + Cf /(Wf C;)]" 2 . 

If the expectation value of this estimator is taken, then g2M,w is replaced with 
(<?2M./i') = /Fu'- This estimator was constructed with the assumption that <?2M 
is scale-invariant, in which case (g2M.1i') = ff^M an d (<72iw) does give the true <?2M- 
However, if g 2 M does vary with scale, then it has to be included inside the integral for 
Df,,, and the estimator's expectation value must be taken as g 2 M at k = fc e ff. We find 
this scale by giving g 2 M(k) a functional form, which we use to find (52A/) = 92M{k c s)- 



+ Although Ref. 6 parametrized the quadrupole asymmetry in terms of g t , not <;2M> this should 
not affect the effective wavenumber of the measurement. 
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We begin with the expectation value of g 2 M,u> , given by 
_ Jg 2 M(k)P ( (k)Q l (k)e ll (k)edk 

{92M ' w} ~ fPcWetQQvW&dk ' ( } 

where P^ is the curvature power spectrum and Qi(k) is the transfer function of the 
CMB temperature fluctuations. Typically, in inflationary models that break scale 
invariance, g 2 M will vary smoothly with In k, and we can write 

92M{k) = Bi + 5 2 ln(fc/fc»), (C.5) 

where fc* is the (arbitrary) pivot wavenumber such that (&*) is constant even when 
the scalar spectral index n s is varied. We then insert this expression into Eq. IC.4l and 
find 

(92M,w) -B 1+ B 2 jp dk)Qi{k)<di , {k)k2dk • (C6) 

The denominator in the second term is just Fw. The numerator can be rewritten by 
using P((k) oc (k/k*) ns such that 

(92M,w)=Bi + B 2 ^^. (C.7) 
on s 

We then insert this expression into the expectation value of Eq. IC.3I and, by using 
(S2m) = Bi + B 2 ln(fc off /fc*), we find 

ln ( fc cff/fc*) = — J7ji ■ (C8) 

We can find the derivative in this expression numerically by using the finite difference 
method with a two-sided derivative centered at the fiducial value n s — 1 with endpoints 
nf = 1 ±0.025. 

We use the fiducial cosmological parameters (except for n s ) to calculate Ci using 
CAMB [53] . along with WMAP's instrumental parameters ot = 1.09 x 10" 5 and 
#fwhm = 21', specifically in the V band. These parameters give us XJs in the range 
of T 2 = 0.99994 to T 400 = 0.17. We select the pivot wavenumber fc» = 0.002 Mpc" 1 . 
This gives us fc° ff MB = 0.0204 Mpc" 1 . 

Finding k c g for the LRG analysis is similar except we include measurements 
of perturbations at 8 redshift slices instead of only one last-scattering surface. By 
making the necessary changes for a galaxy survey analysis, including F l g Z ( ;±2 ) — —C g ; 
(i denotes the redshift slice), we find 

^ lTgilg^g~i MU + j§^{g 2 M,i{i+ 2 ) + ?2Af./(;-2))] ,„ , 

92M - j— ^ — — , (0.9) 

2-n- 2^il,i 11 g,l 

where T g ; = [1 + Cg™/Cg/]~ 2 , C™'/ = Afl/rii, Af2 is the pixel size, and is the 
average number of galaxies per pixel. The expectation value of the estimator g\ M w 
is given by 

^ v _ Jg 2 M(k)P g (k)Wi(k)W{,(k)kidk 

\ 9 ™«'>- / P{k)Wi{k)W},{k)kHk ' 1 j 

where Wf{k) is the (survey) window function. g 2 M(k) again varies smoothly with lnfc, 
which allows us to write 

g 2M (k)=D 1 + D 2 lnk. (C.ll) 
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This parameterization gives us 

/-4 \_ I In kP g (k)Wi(k)Wl,(k)k 2 dk 

{92M,w) ~Di + D 2 f p g{k)W i {k)W i {k)k 2 dk ■ ( C - 12 ) 

The denominator in the second term is just u , , and the numerator can be 
calculated directly. Inserting this into the expectation value of Eq. IC.91 which equals 
L>i + D 2 ln(fc eff ), we find k^ G = 0.151 Mpc" 1 . 
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